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Abstract. Radioactive isotopes such as 44 Ti offer prob- 
ably the most direct probe into nucleosynthesis environ- 
ments in the interior of exploding stars, when the asso- 
ciated 7-ray activities in the explosion remnant are de- 
tected and translated back to the isotopic abundances at 
the time of the explosion. In this paper, we assert that 
the procedure may not necessarily be straightforward, at 
least in the case of 44 Ti, an orbital-electron capture decay 
isotope. Using the analytic model of McKee & Truelove 
(1995) for young supernova remnants, and assuming the 
existence of overdense 56 Fe-dominated clumps that con- 
tain also 44 Ti, we show that a high degree of ionization 
may be caused by the reverse shock so that the electron- 
capture rate of 44 Ti could be significantly reduced from its 
laboratory value. When applied to Cas A, this increases 
under certain conditions the current 44 Ti-activity by a fac- 
tor 1.5 ~ 2.5, which yields a better compatibility between 
the COMPTEL observation of the 1.16 MeV line activity 
associated with the 44 Ti decay and the SN model pre- 
dictions of the initial 44 Ti abundance. This possibility is, 
however, subject to various uncertainties, and in particu- 
lar to the unknown properties and radial distribution of 
the clumps in the ejecta. 

Key words: supernovae: general - supernovae: individ- 
ual: Cas A - ISM: supernovae remnants - Nuclear reac- 
tions, nucleosynthesis, abundances 



1. Introduction 

Exploding massive stars, and supernovae in particular, are 
known to be major sites for the production of a large 
variety of elements heavier than carbon. One of the few 
available ways to study the physics in the deep interior of 
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such stars is the determination of the abundances of sta- 
ble nuclides freshly produced and ejected by the explosion. 
Infrared, optical, and X-ray spectroscopic measurements 
are capable of determining elemental abundances in the 
photosphere during different phases of the outburst. But 
the results from such observations are, in general, sensi- 
tive to the models of line excitations in the photosphere, 
resulting in large correction factors to be applied before 
deducing the isotopic yields at the time of the explosion. 
In addition, results are hampered by the uncertainty of 
the optical depth, and by the possibility of heavy element 
condensation into dust shortly after the explosion, as was 
witnessed in SN 1987A. 

A more direct probe of massive-star interior physics is, 
in principle, to investigate unstable nuclides and to mea- 
sure the 7-rays associated with their (3 decays after they 
have been ejected by the supernova explosion. For an op- 
timum probe, the mean lifetime of such a radioactive iso- 
tope should range from around a few weeks up to about 
10 6 years. The lower limit is set by the requirement that 
the ejecta should become optically thin to 7-rays in a few 
decay times, and the upper limit by instrumental sensitivi- 
ties of 7-ray telescopes (the 7-ray flux from a trace isotope 
must exceed the instrumental noise level, presently in the 
range of 10~ 5 photons cm _2 s _1 , which corresponds, for 
instance, to (several times) 10~ 2 M Q of an intermediate- 
mass isotope with a lifetime of 10 6 y at the Galactic center 
or to 10~ 4 Mq of the isotope in a supernova at a distance 
of a few hundreds of pc) . 

Furthermore, whereas short-lived isotopes will clearly 
trace individual events, long-lived ones with mean life- 
times of the order of 10 6 y or longer will reflect a superposi- 
tion of different supernovae at different times, mixed with 
interstellar matter. Consequently clues to abundances in 
an individual object are only very indirect. 

Only a few isotopes fulfill those constraints (see e.g. 
Diehl and Timmes 1998). Most promising cases are found 
among the Fe group elements, primarily because of their 
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expected large abundances. The 0.847 and 1.238 MeV 7- 
ray lines from the 56 Co — ^ 56 Fe decay (half-life: ti/ 2 = 77 
d) were detected from SN 1987A (Matz et al. 1988; Sandie 
et al. 1988; Mahoney et al. 1988; Rester et al. 1988; Tee- 
garden et al. 1989). There is also evidence for these decay 
lines from the unusually fast and bright Type la supernova 
1991T (Morris et al. 1995, 1997). The 57 Co -^ 57 Fc decay 
(ti/2 = 272 d) is another probe: the 122 and 136 keV lines 
were detected from SN 1987A (Kurfcss et al. 1992; Clayton 
et al. 1992). Cases at the upper end of the favored radioac- 
tive lifetime range are 26 Al (ti/ 2 = 7.4 xlO 5 y) and 60 Fe 
(ii/2 = 1.5 x 10 6 y). The 1.809 MeV line from the 26 Al de- 
cay has been detected and mapped along the entire plane 
of the Galaxy (see review by Prantzos & Dichl 1996). If su- 
pernovae, rather than Wolf Rayet stars, were responsible 
for this 26 Al, the lines from the 60 Fe (3~ decay would be 
expected simultaneously with the 26 Al decay, identifying 
a supernova origin (review by Diehl and Timmes 1998). 
Instrumental sensitivity appears just at the borderline for 
this test. 

In this paper, we focus our discussion on 44 Ti, which 
decays with t x / 2 = 60 y (Ahmad et al. 1998; Gorres et al. 
1998; Norman et al. 1998; Wietfeldt et al. 1999), making it 
ideal for a study of inner-supernova physics within young 
supernova remnants. 44 Ti decays almost uniquely to the 
2 nd excited state of 44 Sc, followed immediately by the al- 
most unique f3 + decay of 44 Sc (ii/ 2 = 4 h) to the 1.156 
MeV excited state of 44 Ca. The 1.156 MeV de-excitation 
line has indeed been observed by the COMPTEL telescope 
on the Compton Observatory from Cas A, a young super- 
nova remnant with an estimated age of 320 y (Iyudin et al. 
1994, 1997; Dupraz et al. 1997). The measured 7-ray flux 
is ~ (4 ± 1) x 10~ 5 photons /cm 2 /s (Iyudin et al. 1997) 
concordant with an upper limit obtained by the OSSE in- 
strument (The et al. 1996). With an adopted distance to 
Cas A of 3.4 kpc (Reed et al. 1995) and the laboratory 
decay rate, the inferred initial mass of 44 Ti is ~ 2 x 10~ 4 
M Q (Iyudin et al. 1997; Woosley & Diehl 1998). 

The current model predictions of the 44 Ti initial mass 
lie in an approximate range of (6 x 10~ 6 ~8x 10 -5 ) M 
for Type-II SNc (Woosley & Weaver 1995; Thielemann et 
al. 1996), and of (3 - 8) x 10" 5 M© for Type-lb SNc 
(Woosley et al. 1995), more or less strongly depending on 
the progenitor masses. Higher values up to 2 x 10~ 4 M Q 
were obtained in some of the Type-II SN models, when 
progenitor masses above 30 M Q combine with high ex- 
plosion energies (Woosley & Weaver 1995), and for a 20 
M Q (SN1987A) model star (Thielemann et al. 1996). Bar- 
ring the possibility that the progenitor of Cas A happened 
to be such a star, one may conclude that COMPTEL ob- 
served significantly more Ti than expected (see, e.g., Fe- 
sen & Becker 1991; Hurford & Fesen 1996 for discussions 
on the progenitor characteristics). 

As far as its /?-decay properties are concerned, 44 Ti is 
a very interesting trace isotope, because its decay mode 
is pure orbital electron capture, which means that fully 



ionized 44 Ti is stable. Even partial ionization of the in- 
nermost electrons should lead to a considerably longer ef- 
fective half-life (Mochizuki 1999). Therefore, the question 
arises whether in supernova remnants 44 Ti could be highly 
ionized and thus more stable for a considerable period of 
time during the evolution. In this case, it would be incor- 
rect to use the half-life measured in the laboratory, and 
initial abundances of 44 Ti as deduced from 7-ray intensi- 
ties could be too high. 

However, there is no simple answer to this question 
and, as we shall emphasize below, the thermodynamic his- 
tory of a remnant has to be known in detail before firm 
predictions can be made. On the other hand, it is inter- 
esting to speculate whether the COMPTEL observations, 
which indicate an amount of 44 Ti in Cas A that appears 
higher than expected, reflect the effects of an increased 
lifetime of the explosively-produced 44 Ti because of tem- 
porary and partial ionization. 

The primary aim of the present paper is to outline pos- 
sible implications of (partial) ionization on the observable 
7-ray line flux from the decay of 44 Ti. We shall present re- 
sults obtained for a variety of different conditions, based 
on a simple model for young supernova remnants which, 
nonetheless, accounts for the features most relevant to this 
question. Of course, our main focus will be on Cas A, the 
best studied case, and the parameters of the model are 
chosen accordingly. 

Section 2 contains an overview of observational aspects 
of SN explosions and remnants which are of relevance in 
the context of this work. In Sect. 3, we describe the em- 
ployed model for young supernova remnants, i.e., the an- 
alytic model of McKee & Truelove (1995), augmented by 
a description of the reverse shock interacting with denser 
cloudlets (Sgro 1973; Miyata 1996). In addition, we de- 
scribe the microphysics used in the model on which the 
calculation of the "effective" decay rate of 44 Ti is based. 
In Sect. 4, we present our results for the time variation of 
the 44 Ti decay rate, the 44 Ti abundance in young super- 
nova remnants and the associated 7-ray activities that can 
be measured. As a specific example, we will preferentially 
use the case of Cas A. Summary and conclusions are given 
in Sect. 5. 

2. Young Supernova Remnants 

Studies of supernovae from the moment of the explosion 
up until they dissolve and merge with the general inter- 
stellar medium has, in general, been guided by observa- 
tional opportunities as combined with model interpreta- 
tions; the interplay of different physical processes, which 
vary spatially and rapidly within a young supernova rem- 
nant, cannot be disentangled from measurements alone. 
Observational windows are (from low to high energy radi- 
ation): Radio maps from electrons synchrotron-emitting in 
magnetic-field structures around the shocked-gas region; 
infrared emission from cool and hot dust within dense 
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clumps embedded in the supernova remnant; optical line 
emission from the edges of dense material embedded in the 
remnant's hot plasma; X-ray emission from the hot, ion- 
ized, gas that has been shocked by the forward blast wave 
and by the reverse shock traveling inward, respectively; 
and 7-rays from long-lived radioactivity such as 44 Ti. 

Here we are mainly interested in the thermal history 
of the radioactive 44 Ti after it has been ejected in a su- 
pernova explosion. 44 Ti as well as other iron-group ele- 
ments are synthesized during the very early moments of 
the explosion of a massive star in the layers adjacent to 
the nascent compact object (a neutron star or black hole). 
Therefore, observations of iron can be used to trace also 
44 Ti. However, not many good cases are known so far. 

SN 1987A is certainly the best studied case. The large 
Dopplcr shifts of iron lines observed in the early spec- 
tra require that iron moves with velocities much higher 
than some of the hydrogen which cannot be explained 
by spherically symmetric explosion models but indicates 
that it is not homogeneously distributed in an expand- 
ing shell but is found in clumpy structures. This con- 
clusion is also supported by the unexpectedly early de- 
tection of X- and 7-rays from the decay of radioactive 
elements (for summaries, see Woosley & Weaver 1994, 
Nomoto et al. 1994), which again suggests that these nu- 
clei have been transported far out into the hydrogen-rich 
shells. Other arguments in favor of this interpretation in- 
clude the smoothness of the light-curve (e.g., Woosley & 
Weaver 1994, Nomoto et al. 1994 and references therein), 
and the time-dependent features of the spectral lines ob- 
served soon after the outburst (Utrobin et al. 1995), e.g., 
in the Bochum event (Hanuschik & Dachs 1987; Phillips 
& Heathcote 1989). In fact multi-dimensional supernova 
simulations (e.g., Herant et al. 1992, 1994; Shimizu et al. 
1994; Shimizu 1995; Burrows et al. 1995; Janka & Miiller 
1995, 1996) have demonstrated that the surroundings of 
the newly-formed neutron star are stirred by hydrody- 
namic instabilities and that inhomogeneities and dumpi- 
ness of the products of explosive nucleosynthesis are likely 
the consequence. However, it cannot be excluded that SN 
1987A is a special case and not comparable with, e.g., 
Cas A. 

44 Ti 7-ray emission is expected from core-collapse su- 
pernovae in general (e.g., Woosley & Weaver 1995; Thiele- 
mann et al. 1996), yet there is a significant deficit in such 
7-ray line sources in the Galaxy for the inferred Galac- 
tic core-collapse supernova rate (Dupraz et al. 1997). The 
44 Ti detection for Cas A (Iyudin et al. 1994; 1997), despite 
the low Fe abundance, therefore has given rise to specula- 
tions about its exceptional nature, particularly because its 
optical and X-ray characteristics support the idea of asym- 
metries and a peculiar circumstellar environment (see be- 
low, and Hartmann et al. 1997 and references therein). It 
was suggested, for example, that the progenitor of Cas A 
might have been a rapidly spinning star, in which case 



more 44 Ti could have been synthesized (Nagataki et al. 
1998). 

Recently, a second Galactic 44 Ti source has been re- 
ported (Iyudin et al. 1998). Its alignment with an also re- 
cently discovered X-ray remnant (Aschenbach 1998) sug- 
gests that it is a very young and most nearby supernova 
remnant, with an age around 680y and a distance of 200pc 
only. This object still provides a puzzle because of the 
absence of radio and optical emission expected for such 
a nearby supernova. Yet, if confirmed, it will provide a 
unique opportunity for the study of supernova-produced 
44 Ti; in which case it is of some interest to see if the mod- 
ified decay rate of ionized 44 Ti addressed in our paper 
could also be important. 

The 44 Ti found in young supernova remnants is prob- 
ably formed during the a-rich freezeout from (near) nu- 
clear statistical equilibrium (Woosley et al. 1973; Woosley 
& Weaver 1995; Woosley et al. 1995; Thielemann et al. 
1996; Timmcs et al. 1996; Timmes & Woosley 1997; The 
et al. 1998). Whereas 44 Ti is obviously fully ionized at 
the time of explosion, 44 Ti ions will recombine with elec- 
trons during the adiabatic cooling phase, which accompa- 
nies the expansion of the exploded star, to become neutral 
after some 1000 s (e.g., Nomoto et al. 1994), a negligible 
timescale when compared with the age of a supernova rem- 
nant. Therefore, during most of this early phase 44 Ti will 
decay with the laboratory rate. The subsequent evolution 
of the remnant may however provide conditions for its re- 
ionization, which is the subject of our modeling effort. 

The evolution of young supernova remnants has been 
modeled extensively, and the spherically-symmetric ex- 
plosion into homogeneous interstellar surroundings ap- 
pears well- understood (McKee and Truelove 1995). A free 
expansion phase is followed by an adiabatic blast-wave 
phase, where interaction with surrounding material pro- 
duces outward and inward moving shock waves leading to 
bright X-ray and radio emission, yet being unimportant 
for the energetics of the remnant. Later phases of signif- 
icant slowing-down and radiative losses of the expanding 
remnant lie beyond the early phase where 44 Ti still decays. 
Young remnants such as Cas A are generally understood 
as being somewhere intermediate between free expansion 
and the second phase, commonly called "Sedov- Taylor" 
phase. 

The evolution of such idealized young supernova rem- 
nants during the ejecta-dominated stage (t < tsr) and the 
Sedov- Taylor phase (t > £st) was described in an analyt- 
ical model by McKee & Truelove (1995). At the interface 
between ejecta and ambient medium a contact discontinu- 
ity occurs, separating the exterior region of shocked and 
swept-up circumstellar gas behind the outward-moving 
blastwave from the supernova remnant interior. Inward 
from the discontinuity, a reverse shock travels through the 
expanding, cold supernova ejecta, heating up the interior 
gas. Bright X-ray emission results from the hot plasma 
on both sides of the contact discontinuity, with higher 
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temperature on the outside heated by the blastwave (> 
keV), as compared with reverse-shock heated gas on the 
inside (< keV; e.g. Vink et al. 1996). Although such a 
two-temperature model appears adequate to describe the 
X-ray emission of many supernova remnants, considerable 
uncertainty remains in detail. For example, the blastwave 
shock rapidly heats the ions entering the blastwave region, 
the thcrmalization times for the X-ray emitting electrons 
may exceed 100 years, so that non-equilibrium models are 
needed for a proper description of X-ray and radio emis- 
sion. The detailed physical conditions within supernova 
remnants are far from being understood, although the 
general evolution follows these fairly simple descriptions 
quite closely, and is controlled by a few parameters, the 
explosion energy, the mass of ejecta, and the surround- 
ing medium density. According to the McKee & Truelove 
(1995) model, an explosion energy in the (1 ~ 3) x 10 51 erg 
range does not seem unreasonable in the case of Cas A for 
the possible ejecta mass of (2 <~ 5) M Q (Tsunemi et al. 
1986; Vink et al. 1996) if the surrounding ambient gas den- 
sity is of the order of 20 cm~ 3 (Tsunemi et al. 1986) and 
the current blastwave radius is (2 <~ 3) pc for an assumed 
distance of 3.4 kpc (Anderson & Rudnick 1995; Holt et al. 
1994; Jansen et al. 1988; Reed et al. 1995). 

Although the general appearance of supernova rem- 
nant images in the radio and in X-rays is that of a 
large-scale shell configuration as expected from the above 
model, additional prominent clumpy structures appear in 
some cases (e.g., Anderson & Rudnick 1995; Koralcsky et 
al. 1998). This suggests that the model outlined so far is 
an oversimplification as far as details are concerned. The 
gross radio and X-ray emissivity may not be very sensi- 
tive to such discrepancy, tracing the electron component 
in the vicinity of the shock region, but the bulk ejected 
mass may be inadequately represented by the inferred 
electron densities and temperatures (e.g., Koralesky et 
al. 1998). Thus, for the Cas A remnant, prominent struc- 
tures have been studied in their respective forms of opti- 
cal knots and filaments, "quasi-stationary flocculi" , "fast- 
moving knots", and "fast-moving flocculi" (e.g., van den 
Bcrgh & Kamper 1983; Reed et al. 1995; Peimbert & van 
den Bergh 1971; Chevalier & Kirshner 1977, 1978, 1979; 
Reynoso et al. 1997; Lagage et al. 1996). There is over- 
whelming evidence for dense structures embedded in ten- 
uous material within the entire remnant, and even outside 
the blastwave shock radius. Obviously the explosion itself 
produces fragments of material, seen now as fast-moving 
knots with their abundance patterns supporting an ejecta 
origin. These clumps might carry heavier elements pref- 
erentially as suggested by observations of fast-moving Fe 
clumps early in SN explosions, e.g. in SN 1987A (e.g., 
Nomoto et al. 1994, Wooden 1997 and references therein), 
but a connection between the instabilities and dumpiness 
early after the supernova explosion and the fragments and 
"bullets" seen in the remnants has not been established 
yet. 



Given all these uncertainties in the evolution of su- 
pernova remnants we do not attempt to model a specific 
object, such as Cas A, in detail here. We rather shall in- 
vestigate by means of an admittedly very simple model, 
varying its parameters within reasonable limits, the poten- 
tial effects of ionization on the 44 Ti abundance estimates 
obtained from 7-ray observations. 

3. The Model 

The remnant model employed in this work uses the an- 
alytic description by McKee & Truelove (1995) for the 
hydrodynamic evolution. The ejecta with mass M e j are 
assumed to be cold and to expand with explosion energy 
E e j into a homogeneous ambient medium which has a hy- 
drogen number density tiro and is composed of 10 hydro- 
gen atoms per helium atom, corresponding to a helium 
mass fraction of 28.6%. Inside the homogeneous ejecta, 
we assume the existence of denser clumps that contain a 
small fraction of the ejecta mass but most of the produced 
44 Ti. The treatment of the interaction of the reverse shock 
with these clumps and of the corresponding effects on the 
44 Ti decay are described below. Ionization of 44 Ti is most 
efficiently induced by electrons. If and how long a high 
degree of ionization is sustained depends on the reverse- 
shock characteristics and the postshock evolution of the 
matter containing Ti. 

3.1. A model for young supernova remnants 

If sufficiently high temperatures happen to be reached be- 
hind the reverse shock, 44 Ti may become fully ionized, in 
which case its decay is prevented until cooling due to the 
expansion of the gas leads to (partial) recombination of 
the innermost electrons. 

To follow the radioactivity of the supernova remnant 
requires knowledge of the evolution of the shocked ejecta. 
The analytical model of McKee & Truelove (1995) yields 
scaled relations for the radius, velocity and postshock tem- 
perature of the reverse shock as functions of time t. Also, 
the model provides the density of the unshocked (homo- 
geneous) ejecta at time t, which allows one to calculate 
the postshock density from the density jump at the re- 
verse shock. In order to describe the density history of 
each mass shell after the reverse shock has passed through 
it, we assume that the shocked matter moves (approxi- 
mately) with the same velocity as the contact discontinu- 
ity. Knowing the motion of the latter, therefore, we can 
estimate the dilution of the shocked gas by the expansion 
of the supernova ejecta. 

The time-dependent analytical solutions are given by 
the three parameters, M e j, E c j and nno- A combination of 
values of these parameters is considered to yield a suitable 
description of a certain supernova remnant, if the model 
gives a present-day radius and velocity of the blastwave 
that is compatible with observational data. 
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We extended the analytic remnant model of McKee & 
Truclovc (1995) by assuming that the iron-group elements 
together with the explosively nucleosynthesized 44 Ti are 
concentrated in overdense clumps of gas. These clumps are 
assumed to move outward through the dominant mass of 
homogeneous supernova ejecta with high velocity. Since 
the clumps should contain only a minor fraction of the 
total mass of the ejecta, the dynamics of the young su- 
pernova remnant as described by the McKee & Truelove 
(1995) model is assumed not to be altered by their pres- 
ence. The radial position r of the iron clumps within the 
ejecta can be measured by the relative mass coordinate q 
which is defined as 



47rr 3 p cj 
3M ei 



(1) 



3.2. Thermodynamic conditions just behind the reverse 
shock 

We now consider the properties of a clump of metal-rich 
matter within the homogeneous ejecta. Denoting the den- 
sity enhancement factor of the assumed clump relative to 
the surrounding ejecta by a c impj we obtain the clump den- 
sity after the reverse shock has passed at time t s h as 



p(t s h) = 4a c i m pp j(i sh ) , 



(2) 



for an adiabatic index 7 = 5/3. 

Correspondingly, the post-reverse shock temperature 
of the clump is related to that of the surrounding uniform 
ejecta, T cj (i sh ), as 



T(t sh ) 



/3(a c im P ) M*sh) 



C^clmp Mej 



T j(t s h) , 



(3) 



which generalizes by the factor /z(i s h)//x e j the result of 
Sgro (1975) to the case considered here, where the chem- 
ical composition of the clump is different from that of 
the surrounding ejecta as described by the correspond- 
ing mean molecular weights /u(i s h) and /i c j, respectively. 
r e j(i sh ) is given by 



T cj (U) = ^ r 2 (U), 



(4) 



where v r is the reverse-shock velocity in the rest frame 
of the unshocked ejecta (McKee & Truelove 1995). In 
Eq. (3), (3 is the ratio of the pressure values in the re- 
gions behind the reflected shock and in front of it. The 
reflected shock is formed when the reverse shock hits a 
dense clump in the ejecta (Sgro 1975; Miyata 1996). The 
ratio (3 is related to a c i mp as 



^clmp 



1 



l-(3 



1) 



(5) 



for 7 = 5/3 (Sgro 1975). In particular, f3 = 1 if a c imp = 1 
(no dumpiness), and /3/<*cimp appearing in Eq. (3) is less 
than unity for a c i mp > 1. 

3.3. Post reverse-shock evolution 

We now follow the fate of the homogeneous ejecta and 
of the matter in clumps after they have experienced the 
impact of the the reverse shock at time t s h- As we have 
mentioned earlier, we assume that the shock itself would 
continue to proceed inward through the cold ejecta as pre- 
scribed by the McKee & Truelove (1995) model. The ther- 
modynamic and chemical evolution of the clumps at times 
t > tgh is treated as described below. 

Since we trace the clump material, it is preferable to 
introduce Lagrangian particle abundances, defined by 



where p c j is the uniform background density of the ejecta 
which expand homogeneously. 71 = 



dN 



(6) 



where dN is the number of particles at the mass coordi- 
nate q in the interval dq. 

3.3.1. Ionic abundances 

At a given time £(> t s h), the number abundance m of a 
nuclear species k with charge Z and i orbital electrons (i.e. 
Z — i times ionized, and rik = ^2 i>0 ni) follows 



dn t 



A, 



ion, i 

+^ion,i+l^i+l + A 



rcc.z— l^i— 1 •> 



(7) 



where Xpjfy), \ion,i(t), and A reCi j(i) arc the nuclear f3- 
decay, ionization and recombination rates, respectively, of 
the ion. For 44 Ti, \p ^ 0, and thus n(t > i s h) < n(t s h) = 
n(0)exp(— Ai a btsh), where A; a b is the 44 Ti decay rate in the 
cold ejecta, i.e. its laboratory value. 

Since we assume that material which contains 44 Ti 
consists dominantly of one single stable nuclide ( 56 Fe in 
practice), the number of ionization electrons, n c , is given 
by ^2iy (Z — i)rii for that species. In practice, we further 
assume that both 56 Fe and 44 Ti are ionized to some degree 
in the reverse shock itself, which we cannot resolve. Typi- 
cally, we start our network calculations with ionic systems 
of 10 bound electrons, This choice is made since our ma- 
jor concern is to see the reduction of the orbital-electron 
capture rates in the domains of high shock temperatures. 

The reaction rates entering in Eq. (7) are discussed in 
Appendix A. 

3.3.2. Ion and electron temperatures 

The network requires the knowledge of the thermody- 
namic evolution of the shocked matter, and in particular 
of the electron temperature, T e , which defines the Boltz- 
mann distribution, and of the matter density p (or the 
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ion number density, n; n), which enters in the absolute 
electron number density n e . 

The decrease by expansion of the temperature T(t) 
from T(t s h) of Eq. (3) would roughly be proportional to 
[p{t)j p{t s h)] 2 ^ 3 (for 7 = 3/5 and t > t sh ). Two considera- 
tions are due, however. One concerns the possibility that 
both T c and the ion temperature, T ion , may not necessar- 
ily be the same as T. Secondly, some non-adiabatic effects 
may appear as the result of ionization and recombination. 

It is sometimes asserted that the collisionless equilibra- 
tion between Tj on and T c would occur quickly, in which 
case Tion = T c = T already at t = t sn . Save this possi- 
bility, the cquipartition of the shock energy to ions and 
electrons may not be achieved at once. The equilibration 
process through the Coulomb interaction is described by 
a non-linear equation (Spitzer 1962), 



die 
dt 



(8) 



in terms of the equilibration time r eq (t), which is a func- 
tion primarily of Tj on , T c and the number densities of ions 
and electrons. Assuming that T c <C T lon at t = t sn , and 
that the total kinetic energy 



Rb and velocity iv Below we consider mainly the parame- 
ter space which is consistent with these constraints. While 
we treat a c imp as a free parameter, a value of the order of 
10 seems to be reasonable as a first guess, based on the 
density contrast which develops in hydrodynamic insta- 
bilities in the envelope of the progenitor star during the 
supernova explosion (Fryxell et al. 1991). 
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is constant during a short time interval St, we solve Eq. (8) 
simultaneously with the abundance network equations to 
determine the degree of ionization of a dominant species, 
and thus n c . 

In order to include the non-adiabatic effect caused by 
ionization and recombination, we assume that the energy 
required for unbinding additional electrons (net energy 
loss) is taken from electrons only. Namely, after perform- 
ing the network calculation at each time step and solving 
Eq. (8), we make a replacement 



T e exp 



2 SQ 



Sn e _ 
n e 3 n G T G * 



(10) 



where 5n e and SQ are the change of n c during St and the 
corresponding net energy consumption. The equilibrium 
temperature T is then calculated from Eq. (9). 

Finally, T(t+St), T ion {t+St), and T c {t+St) are obtained 
from the current values by multiplying [p(t + St)/ p(t)] 2 ^ 3 . 

4. Results 

We present here our results, which we obtained with the 
model described above. They depend on the following pa- 
rameters: E e j , M e j , tiho 7 Q^cimp an d two values for q, given 
as q and qi, which define the lower and upper bound- 
aries of the mass shell where the 44 Ti-containing clumps 
are assumed to exist. As we mentioned earlier, the values 
of the first three can be constrained in the case of Cas A 
from the observational analyses of the blast-wave radius 



Fig. 1. Examples for the evolution with time t (in years) of the 
56 Fe number density riFe {top), and of the electron temperature 
T e (bottom) in the 56 Fe-dominated clumps at four different lo- 
cations indicated by the mass coordinate q. The dashed lines 
represent the asymptotic behavior of the equilibrium tempera- 
tures T(t) (as described in Sect. 3.3.2) at q = 0.4 and 0.6. The 
left edge of each solid curve corresponds to t = i s h- For exam- 
ple, the reverse shock took 50 y to reach q = 0.6. The results are 
given for the following set of parameter values: E c j = 3x 10 51 
erg, M e j = 3 Mq, nm = 15 cm~ 3 and a c im P = 10 



4-1. Thermodynamic evolution of post-shock clumps 

Figure 1 shows the evolution of the 56 Fe number den- 
sity riFe and of the electron temperature T e in the 56 Fe- 
dominatcd clumps at four different locations in the ejecta, 
given by the corresponding mass coordinate q, after the 
reverse shock has passed the clumps. The post-shock ex- 
pansion decreases the density as well as the equilibrium 
temperature, while T c first increases in its attempt to 
reach equilibrium between electrons and ions. The results 
are plotted for the chosen set of parameter values quoted 
therewith. 

4-2. Effective decay rates of 44 Ti and its radioactivity as 
an observable 

Let us now consider 44 Ti in a clump that is located any- 
where in the mass coordinate < q < 1 of the ejecta. If 
44 Ti is near the surface (large q), it suffers from the re- 
verse shock early. Since the electron temperature is still 
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Fig. 2. Top panel: The average numbers N e of electrons that 
arc bound to 44 Ti and to 56 Fe at mass coordinate q and at time 
t (in y) after the explosion. Middle: Effective decay rates A post sh 
of the shocked 44 Ti during the time span between i sh and t 
[Eq. (11)] in units of the laboratory rate Ai a b- Bottom: The 
corresponding 44 Ti radioactivity observable at time t ("age") 
relative to the case that assumes no reduction of the /3-decay 
rate. The parameter values used here are the same as those for 
Fig. 1 



relatively low (see Fig. 1), there is no possibility to reach 
a high degree of ionization (thus an appreciable reduction 
of the decay rate). On the other hand, 44 Ti embedded in 
clumps located in the innermost region may become highly 
ionized as the temperature is generally highest, resulting 
in the smallest 44 Ti decay rate. Since i s h is closer to t, 
the "age" of the remnant, however, the net effect of the 
reduced decay rate appears weaker. Therefore, the maxi- 
mum effect of the retarded /3-decay is obtained if 44 Ti is 
at intermediate values of q. This situation is depicted in 
the panels of Fig. 2. 

The top panel of Fig. 2 shows the average numbers 
of electrons that are bound to 44 Ti and to 56 Fe at t = 
100, 200 and 300 years after the explosion. The increase of 
N e toward the low side of q indicates that the time span 
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Fig. 3. 44 Ti radioactivities at times t (in y) in the clumps 
located at radial positions q, relative to the radioactivity 
of the model with most detailed input physics, when the 
non-adiabaticity of the thermal evolution of the clumps is ig- 
nored (upper panel) and when instantaneous equilibration be- 
tween T e and T lon (toward T) is assumed (lower panel) 



between i s h and t was not long enough for ionization. At 
the high end of q- values, the effects of post-shock recombi- 
nation are visible. When combined with n-F e from Fig. 1, 
iV e ( 56 Fe) determines the total number of ionization elec- 
trons, n e (t), in a clump: n c = [Z — AT c ( 56 Fe)] np c with 
Z = 26. 

The retardation of the 44 Ti decay by ionization is illus- 
trated in the middle panel of Fig. 2. The effective decay 
rate in the post-shock period, A posts h, is defined as the 
time-average of A c ff [given by Eq. (18) in Appendix A] be- 
tween the time i s h and a time t > i s h- It can be expressed 
as 



A 



1 



postsh 



t — ish 



In 



n(t) 



n(t sh ) 



(11) 



where n(t) is the 44 Ti abundance at mass coordinate q and 
time (age) t, and n(t s h) is the corresponding value at the 
shock impact time t s h, which has been reduced from its 
initial value n(0) according to n(t s h) = n(0)exp(— Ai a bi s h)- 
In order to estimate the observable consequences of 
the reduced /3-decay rates, one has to recall that the mea- 
surable 7-ray activity per (normalized) unit mass of the 
remnant is the product of the current /3-decay rate and 
the current abundance, n(t)X e s(t). This introduces a non- 
linear relation between the observable activity given in 
the bottom panel of Fig. 2 and the effective decay rates 
displayed above, particularly in the absence of quick re- 
combination. 

In Fig. 3, we compare the results obtained with the 
input physics as described in Sect. 3 to the cases where 
either the non-adiabaticity of the equilibration process of 
electrons and ions owing to the energy consumption by 
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the ionization process is neglected (upper panel), or where 
instantaneous equilibration (i.e., T = Ti on = T always) is 
assumed (lower panel). It can be seen that in both cases 
the trends of the "full" model shown in Fig. 2, where both 
the non-adiabaticity and the gradual equilibration process 
are taken into account, are somewhat enhanced. This can 
be understood by the fact that the omission of cither of 
the two effects leads to higher electron temperatures and 
thus to a higher degree of ionization in the clumps. 

4-3. 44 Ti radioactivity in Cas A 

The 44 Ti radioactivity of the whole supernova remnant as 
observable at an age t, relative to the activity level that 
would result in the absence of any retardation of the /3- 
decay rate by ionization, is defined by the ratio 



91 



90 



n(q,t)X eS (q,t)dq / iV(0)exp(- W)Ai ab . (12) 



Here, N(0) is the total number of 44 Ti nuclei in the 56 Fe 
clumps at t = 0. For a remnant like Cas A, which is in the 
transition phase from the ejecta-dominated stage to the 
Scdov- Taylor phase, the reverse shock has passed through 
most of the ejecta. Therefore, the factor A must be ex- 
pected to be larger than unity when the remnant is ob- 
served at this stage of its evolution. This implies that the 
initial abundance inferred from the current 7-ray activity 
due to the 44 Ti decay is lower by a factor A^ 1 than the 
amount of produced 44 Ti that is estimated on the basis of 
the laboratory decay rate Ai a b- 

For certain combinations of the characterizing param- 
eters of the supernova remnant model, E e - } , M e j and nno; 
the delayed decay of ionized 44 Ti can have a sizable effect 
on the 7-ray activity of Cas A which is observable at the 
present age of about 320 years. In Tables 1-3 the ratio 
A denotes the currently observable activity according to 
our model, normalized to the expected activity if ioniza- 
tion effects are not taken into account. The results depend 
on the assumed density enhancement a c i mp in the clumps 
relative to the density of the homogeneous ejecta, and on 
the location of the clumps within a shell in the expanding 
ejecta bounded by the lower and upper mass coordinates 
qo and qi, respectively. In addition, in Tables 1-3 numbers 
are given for the blast-wave radius R\, and the blast-wave 
velocity Vh as predicted by the McKee & Truelove (1995) 
model at the current age of the Cas A supernova remnant. 

The ratio A exhibits the following tendencies. For fixed 
explosion energy and ejecta mass, A increases with the 
density of the circumstellar gas of the supernova remnant 
because the expansion and dilution of the ejecta are slowed 
down for higher nno (he., t>b and R\ are smaller at the 
same time). An increase of the ratio A can also be seen 
when the explosion energy becomes larger but all other 
parameters are kept constant. The opposite trend is visible 
if the explosion energy and ambient density are fixed, in 
which case the ratio A decreases with larger ejecta mass. 



Table 1. Illustration of the variability of the 44 Ti radioac- 
tivity ratios A as defined by Eq. (12) with different values of 
the model parameters. The tabulated results are at the time 
(age) of 320 y after the explosion for an explosion energy of 
E c j = 1 x 10 51 erg and various combinations of values for the 
ejecta mass M c j , the hydrogen number density in the ambient 
medium, nno, the over-density factor of the clumps, a c im P , and 
the lower and the upper boundaries qo and q\ of the mass shell 
where the clumps are assumed to exist. Rb and Vb are the ra- 
dius and velocity of the blast-wave at that time as computed 
from the McKee & Truelove (1995) model for the given sets of 
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Table 2. Same as Table 1 but for E cj = 2 x 10 jl erg Table 3. Same as Table 1 but for £ c j = 3 x 10 jl erg 





nHo 


Rb 


Vb 




go 


<?i 


A 


Mej 


nHo 


Rb 




Q^clmp 


go 


gi 


A 


(M ) 


(cm" 3 ) 


(PC) 


(10 3 km/s) 










(M ) 


(cm- 3 ) 


(PC) 


(10 3 km/s) 










2 


5 


2.44 


3.63 


5 


0.4 


0.6 


1.19 


2 


5 


2.68 


3.84 


5 


0.4 


0.6 


1.22 


2 


5 


2.44 


3.63 


10 


0.0 


0.5 


1.03 


2 


5 


2.68 


3.84 


10 


0.0 


0.5 


1.04 


2 


5 


2.44 


3.63 


10 


0.4 


0.6 


1.33 


2 


5 


2.68 


3.84 


10 


0.4 


0.6 


1.36 


2 


5 


2.44 


3.63 


10 


0.2 


0.8 


1.26 


2 


5 


2.68 


3.84 


10 


0.2 


0.8 


1.36 


2 


5 


2.44 


3.63 


20 


0.4 


0.6 


1.51 


2 


5 


2.68 


3.84 


20 


0.4 


0.6 


1.70 


2 


15 


2.01 


2.80 


5 


0.4 


0.6 


1.53 


2 


15 


2.20 


3.00 


5 


0.4 


0.6 


1.44 


2 


15 


2.01 


2.80 


10 


0.0 


0.5 


1.14 


2 


15 


2.20 


3.00 


10 


0.0 


0.5 


1.09 


2 


15 


2.01 


2.80 


10 


0.4 


0.6 


2.07 


2 


15 


2.20 


3.00 


10 


0.4 


0.6 


2.13 


2 


15 


2.01 


2.80 


10 


0.2 


0.8 


1.53 


2 


15 


2.20 


3.00 


10 


0.2 


0.8 


1.66 


2 


15 


2.01 


2.80 


20 


0.4 


0.6 


2.03 


2 


15 


2.20 


3.00 


20 


0.4 


0.6 


2.57 


2 


30 


1.77 


2.40 


5 


0.4 


0.6 


2.12 


2 


30 


1.93 


2.57 


5 


0.4 


0.6 


1.87 


2 


30 


1.77 


2.40 


10 


0.0 


0.5 


1.36 


2 


30 


1.93 


2.57 


10 


0.0 


0.5 


1.22 


2 


30 


1.77 


2.40 


10 


0.4 


0.6 


2.67 


2 


30 


1.93 


2.57 


10 


0.4 


0.6 


2.85 


2 


30 


1.77 


2.40 


10 


0.2 


0.8 


1.76 


2 


30 


1.93 


2.57 


10 


0.2 


0.8 


1.87 


2 


30 


1.77 


2.40 


20 


0.4 


0.6 


2.11 


2 


30 


1.93 


2.57 


20 


0.4 


0.6 


2.77 


3 


15 


1.96 


2.90 


5 


0.4 


0.6 


1.47 


3 


15 


2.16 


3.08 


5 


0.4 


0.6 


1.47 


3 


15 


1.96 


2.90 


10 


0.0 


0.5 


1.09 


3 


15 


2.16 


3.08 


10 


0.0 


0.5 


1.09 


3 


15 


1.96 


2.90 


10 


0.4 


0.6 


1.74 


3 


15 


2.16 


3.08 


10 


0.4 


0.6 


2.01 


3 


15 


1.96 


2.90 


10 


0.2 


0.8 


1.33 


3 


15 


2.16 


3.08 


10 


0.2 


0.8 


1.49 


3 


15 


1.96 


2.90 


20 


0.4 


0.6 


1.54 


3 


15 


2.16 


3.08 


20 


0.4 


0.6 


2.01 


3 


30 


1.74 


2.46 


5 


0.4 


0.6 


2.04 


3 


30 


1.90 


2.63 


5 


0.4 


0.6 


2.07 


3 


30 


1.74 


2.46 


10 


0.0 


0.5 


1.30 


3 


30 


1.90 


2.63 


10 


0.0 


0.5 


1.30 


3 


30 


1.74 


2.46 


10 


0.4 


0.6 


2.13 


3 


30 


1.90 


2.63 


10 


0.4 


0.6 


2.63 


3 


30 


1.74 


2.46 


10 


0.2 


0.8 


1.53 


3 


30 


1.90 


2.63 


10 


0.2 


0.8 


1.72 


3 


30 


1.74 


2.46 


20 


0.4 


0.6 


1.57 


3 


30 


1.90 


2.63 


20 


0.4 


0.6 


2.12 


4 


15 


1.92 


3.00 


5 


0.4 


0.6 


1.36 


4 


15 


2.12 


3.16 


5 


0.4 


0.6 


1.42 


4 


15 


1.92 


3.00 


10 


0.0 


0.5 


1.04 


4 


15 


2.12 


3.16 


10 


0.0 


0.5 


1.06 


4 


15 


1.92 


3.00 


10 


0.4 


0.6 


1.48 


4 


15 


2.12 


3.16 


10 


0.4 


0.6 


1.78 


4 


15 


1.92 


3.00 


10 


0.2 


0.8 


1.20 


4 


15 


2.12 


3.16 


10 


0.2 


0.8 


1.34 


4 


15 


1.92 


3.00 


20 


0.4 


0.6 


1.31 


4 


15 


2.12 


3.16 


20 


0.4 


0.6 


1.64 


4 


30 


1.71 


2.53 


5 


0.4 


0.6 


1.82 


4 


30 


1.88 


2.68 


5 


0.4 


0.6 


2.02 


4 


30 


1.71 


2.53 


10 


0.0 


0.5 


1.20 


4 


30 


1.88 


2.68 


10 


0.0 


0.5 


1.26 


4 


30 


1.71 


2.53 


10 


0.4 


0.6 


1.74 


4 


30 


1.88 


2.68 


10 


0.4 


0.6 


2.26 


4 


30 


1.71 


2.53 


10 


0.2 


0.8 


1.36 


4 


30 


1.88 


2.68 


10 


0.2 


0.8 


1.55 


4 


30 


1.71 


2.53 


20 


0.4 


0.6 


1.35 


4 


30 


1.88 


2.68 


20 


0.4 


0.6 


1.71 


5 


15 


1.87 


3.10 


5 


0.4 


0.6 


1.25 


5 


15 


2.08 


3.24 


5 


0.4 


0.6 


1.35 


5 


15 


1.87 


3.10 


10 


0.0 


0.5 


1.01 


5 


15 


2.08 


3.24 


10 


0.0 


0.5 


1.02 


5 


15 


1.87 


3.10 


10 


0.4 


0.6 


1.31 


5 


15 


2.08 


3.24 


10 


0.4 


0.6 


1.57 


5 


15 


1.87 


3.10 


10 


0.2 


0.8 


1.12 


5 


15 


2.08 


3.24 


10 


0.2 


0.8 


1.24 


5 


15 


1.87 


3.10 


20 


0.4 


0.6 


1.20 


5 


15 


2.08 


3.24 


20 


0.4 


0.6 


1.42 


5 


30 


1.68 


2.59 


5 


0.4 


0.6 


1.61 


5 


30 


1.85 


2.74 


5 


0.4 


0.6 


1.88 


5 


30 


1.68 


2.59 


10 


0.0 


0.5 


1.11 


5 


30 


1.85 


2.74 


10 


0.0 


0.5 


1.19 


5 


30 


1.68 


2.59 


10 


0.4 


0.6 


1.49 


5 


30 


1.85 


2.74 


10 


0.4 


0.6 


1.94 


5 


30 


1.68 


2.59 


10 


0.2 


0.8 


1.23 


5 


30 


1.85 


2.74 


10 


0.2 


0.8 


1.41 


5 


30 


1.68 


2.59 


20 


0.4 


0.6 


1.24 


5 


30 


1.85 


2.74 


20 


0.4 


0.6 


1.47 



10 



Y. Mochizuki et al.: Ti: Its effective decay rate in young supernova remnants, and its abundance in Cas A 



Although the blast-wave velocity at the present time is 
higher, it was lower early after the explosion and therefore 
the blast-wave radius is smaller for the larger ejecta mass. 

It is not easy to give simple explanations for these 
trends. The ratio A, which scales with the current abun- 
dance of 44 Ti and its effective decay rate, reflects the whole 
time evolution of the remnant and contains contributions 
from all parts of the shocked ejecta with their different ion- 
ization histories. The nonlinear dependence of A on the 
remnant parameters is a result of the interplay between 
a number of effects. For example, the post-shock density 
and temperature in the ejecta, in particular in the clumps, 
are important for the degree of ionization. The time when 
the reverse shock hits the clumps determines the time left 
for the equilibration between electrons and ions and for 
the duration of the delay of the 44 Ti decay. This delay 
lasts from the moment of nearly complete ionization until 
recombination takes place again. 

In order to obtain a large effect on the 44 Ti radioac- 
tivity that is measured in Cas A presently, a high post- 
shock density and temperature are favorable for efficient 
ionization of 44 Ti. On the other hand, the present-day 7- 
ray emission of Cas A means that the temperature must 
not have decreased too slowly and the density not too 
rapidly. If that happened, 44 Ti would be hindered from 
recombination and the observed 7-ray emission could not 
be explained. The maximum effect from the delayed de- 
cay of ionized 44 Ti is found if most of the 44 Ti-carrying 
clumps have been mixed roughly half-way into the ejecta, 
i.e., if they are assumed to be located in the q-intcrval 
between 0.4 and 0.6 (compare Fig. 2), and if the density 
enhancement in the clumps is around 10 for the higher 
explosion energies of E C j — 2 x 10 51 erg and 3 x 10 51 erg 
(Tables 2 and 3, respectively). In case of the lowest con- 
sidered value of the explosion energy, E C j = 1 x 10 51 erg 
(Table 1), a clump over-density of about 5 seems prefer- 
able because of the higher density of the more slowly ex- 
panding ejecta. Higher than optimum a c i mp -values do not 
ensure rapid and strong ionization, whereas smaller than 
optimum a c i mp -values do not allow for quick recombina- 
tion so that the decay rate stays low. If 44 Ti is located 
very far out in the remnant (large q), the temperature be- 
hind the reverse shock is too low for efficient ionization 
(see Fig. 1). If 44 Ti sits deep inside the ejecta (small q), 
the clumps have been reached by the reverse shock not 
sufficiently long ago so as to show a significant effect from 
44 Ti ionization. 

We found maximum values of A between roughly 1.5 
and 2.5 for a large number of different combinations of 
remnant parameters and assumed clump locations and 
density enhancement factors. This means that the present- 
day 44 Ti radioactivity of Cas A could be 50% up to more 
than a factor of 2 higher for a certain amount of Ti, if 
the decay half-life of 44 Ti was stretched by ionization ef- 
fects. Or, reversely, the current 7-ray emission of Cas A 
due to 44 Ti decay might be explained with a significantly 



lower production of this nucleus during the supernova ex- 
plosion. Figure 4 shows that the relative effect from the 
retardation of the 44 Ti decay will increase with the age of 
the remnant. 

For the different considered combinations of parame- 
ter values, Cas A has reached slightly different stages of 
the remnant evolution around the transition time from 
the ejecta-dominated to the Sedov- Taylor phase. The tab- 
ulated blast-wave radius and velocity can be compared 
with observations which give values for the circumstellar 
density of about 20 hydrogen atoms per cm 3 , an estimated 
ejecta mass of about (2 <~ 5) M Q and a blast-wave radius 
of 2-3 pc at a distance of 3.4 kpc (see references given in 
Sect. 2). 

5. Summary 

We examined the effects of the reverse shock on the 44 Ti 
/3-decay rate in young supernova remnants. For this pur- 
pose, we employed the analytic remnant model of McKee 
& Truelove (1995), which was used to describe the hydro- 
dynamic evolution. We assumed that 44 Ti is carried by 
56 Fe-dominatcd, dense clumps which are mixed into the 
otherwise homogeneous supernova ejecta and account for 
a minor fraction of the ejecta mass. Strong ionization of 
Ti due to the heating by the reverse shock and a siz- 
able influence on the 44 Ti decay is obtained if the clumps 
are assumed to be located at intermediate values of the 
radial mass coordinate q. The observationally important 
change of the 44 Ti decay activity, however, depends also 
on the values of the remnant parameters, namely the ex- 
plosion energy, the ejecta mass, the density of the ambient 
medium and the assumed overdensity factor of the clumps 
relative to the embedding homogeneous ejecta. 

We applied our model to the case of Cas A, a young 
supernova remnant with an age of about 320 years. Ob- 
servational data are considered to constrain the parame- 
ter space which was explored. We found that under cer- 
tain conditions the ionization of 44 Ti and the correspond- 
ing delay of its decay can yield an up to three times 
higher 44 Ti activity at the present time than predicted 
on the grounds of the laboratory decay rate. This effect 
is large enough to reduce the apparent discrepancies be- 
tween the 44 Ti production in the explosion as inferred from 
the COMPTEL 7-line measurements and the theoretical 
expectations from the current supernova nucleosynthesis 
models. 

We emphasize again that to get an enhanced 44 Ti de- 
cay activity we had to assume 44 Ti to be associated with 
inhomogeneities in the ejecta which contain the Fe-group 
elements. Hydrodynamical models have so far not been 
able to predict the fraction of Fe and 44 Ti in such clumpy 
structures, the overdensity of the clumps relative to the 
surrounding material at the considered age of a super- 
nova remnant, the dynamics of the clumps, and their size 
and distribution. Multi-dimensional supernova models are 
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Fig. 4. The decrease with time of the amount of 44 Ti in the 
supernova remnant, N(t), by /3-decays, starting from the ini- 
tial amount iV(0) (right scale in logarithm), and the corre- 
sponding, mostly increasing activity ratio A (left scale). Val- 
ues of 3 M (upper panel) and 5 M Q (lower) are used for M e j, 



erg, «ho = 15 cm and a c in 
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adopted in both cases. The different lines show results obtained 
for different assumptions of the radial locations of the clumps: 
q — [qo,qi] = [0.0, 0.5] (solid line), [0.4, 0.6] (dash-dotted line), 
and [0.2, 0.8] (dashed line). The exponential decrease of the 
44 Ti abundance with the laboratory rate is plotted by the dot- 
ted line 



ties concerning the nuclear reaction rates which are im- 
portant to determine the 44 Ti yield (The et al. 1998). 

On the other hand, more detailed X-ray and 7-ray ob- 
servations with a better spatial (and energy) resolution are 
desirable to confirm or reject the implications of the de- 
scribed model, e.g., the existence of Fe-containing knots 
with temperatures in the approximate range of (1 <~ 5) 
keV (cf. Tsunemi 1997). The current angular resolution of 
the ASC A/SIS observations, about one arcmin, is insuf- 
ficient to measure the thermodynamic properties of indi- 
vidual clumps in the ejecta of Cas A. Only information is 
available for the average electron temperature in the mat- 
ter that is heated by the forward and reverse shocks. Com- 
paring the observations with predictions from the remnant 
model used in this paper would therefore require moni- 
toring of the thermodynamic history of the homogeneous 
ejecta. For this, one would need to compute the thermal 
state of the electrons, and would have to specify the un- 
known composition of the ejecta, which is very uncertain 
because it depends on the type of progenitor and explo- 
sion. The described investigations of 44 Ti decay in super- 
nova remnants could be done without such additional as- 
sumptions, but the models do not provide a description 
of the thermodynamic state of the remnant, which allows 
for a direct comparison with current observational data. 

Finally, we mention that we applied our analysis also to 
the new supernova remnant, RX J0852. 0-4622, discovered 
in the direction of the Vela remnant by Aschenbach (1998), 
in which a source of 44 Ti line emission has been detected 
(Iyudin et al. 1998). However, for the estimated low den- 
sity of the ambient medium of less than 0.04(c?/500pc) -0 ' 5 
cm~ 3 with a distance d to the remnant (Aschenbach 1998), 
we found that the reverse shock does not heat the ejecta 
to sufficiently high temperatures to ionize Ti. 
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Appendix A: reaction rates and related uncertain- 
ties 



called for, which connect the very early phase of the ex- 
plosion with the remnant evolution a few hundred years 
later. Self-consistent multi-dimensional simulations have 
yet to be carried out for an improved theoretical picture 
of the explosive 44 Ti production in supernovae. All pub- 
lished theoretical yields have so far only been obtained 
with ad-hoc assumptions, most often by spherically sym- 
metric models and/or simulations of the supernova explo- 
sions which were started from artificial initial conditions 
("piston" models) rather than self-consistent situations. 
We add that such calculations also suffer from uncertain- 



We estimate here the /3-decay, ionization and recombina- 
tion rates, which appear in Eq. (7) . We also evaluate their 
uncertainties in relation to the results presented in Sect. 4. 

The orbital- electron capture rate of 44 Ti: We first for- 
mulate the /3-decay rate of (near) neutral 44 Ti. The cor- 
responding half- life t 1 / 2 has been well determined to be 
60 y (Ahmad et al. 1998; Gorres et al. 1998; Norman et 
al. 1998; Wietfeldt et al. 1999). The 44 Ti /3-decay is by 
orbital-electron capture almost exclusively to the 146 keV 
excited (0~) state of 44 Sc (Endt 1998), the corresponding 
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Q-value being 266 keV. We thus write the laboratory rate, 
Ai a b = In 2/t 1/2 , as 



Alab = Ak + Ali + .... 



(13) 



where Ak, Ali etc. refer to the rates for partial captures of 
K (lsi), LI (2s i) electrons etc. In the "normal" approx- 
imation (e.g. Konopinski & Rose 1965), Ak for example 
reads 



9l 

4tt 2 



\Mo\ 2 (Q + B K ) 2 g 2 K 



(14) 



for the transition of our current interest, where Ga is the 
axial- vector coupling constant of weak interaction, Mq is a 
linear combination of two major nuclear matrix elements 
in the rank-0 first-forbidden (0 + — ► 0~) transition, Q is 
the Q-value, and Bk is the binding energy of the K-shell. 
Thus, q v = Q + Bk is the energy available to the emitted 
neutrino. (The nuclear recoil energy is very small and is 
ignored here.) Finally, gn is the larger component of the 
radial wave- functions of the K-electron, which is evaluated 
at the nuclear radius. 

Since Q is much larger than the binding energies 
(about 5 keV for K-electrons), the relative ratios of the 
partial decay rates become nearly equal to those of the re- 
spective wave-functions squared. Considering in addition 
that aZ is not very large, we have 



Ali ~ o Ak 



(15) 



since gf^ ~ 9k m the \ow-aZ approximation. With the 
use of the same approximation, one finds that the capture 
rates of the higher-shell electrons are small. 

Correspondingly, the orbital-electron capture rates of 
44 Ti with i orbital electrons [i.e., (22-i)-times ionized] arc 
approximately given by 



A/3,i ~ Pi, kK,K +Pi,LlAi,LI 



(16) 



where p^K and p^li are the occupancies ([0,1]) of the K 
and Li orbits, respectively, and A^k and A^li are the cor- 
responding capture rates for the filled shells. We approxi- 
mate the partial capture rates by 



A, 



-A, 



K+LI 



-A 



lab ■ 



A, 



1 



1 



LI 



Aj,K+LI ~ 7rAi a b , (17) 



9 ' " 9 ' 9 ' ^ 9 

so that the total (effective) decay rate reads 



A e ff 



8 17 
9 n2 + 18" 3 



E 

i>4 



n^Aiab / y"/ 



(18) 



with those approximations, including the neglect of the 
decrease of screening and the slightly altered energetics 
along with ionization, are small enough (less than 10 %) 
to be ignored in the present work. 

The ionization rates of 44 Ti and 56 Fe: The rates of 
electron-induced ionization are given by 



A;, 



E 

i>l 



(Pi* ( a ion,i V ) + ^> LI (°"io L n,V 



,(19) 



where n e is the electron number density and (ai on ,iv) is the 
Maxwellian- average of the cross section (times the electron 
velocity v) for ionization of filled K-, LI-... shell electrons. 
As in Eq. (16), p^K and p^Li stand for occupancies of the 
K- and Li shells in the i-th ion. Our estimates of 0i on ,i 
rely on the experimental data on K-shell ionization by 
electrons in the relevant energy range of 6 — 50 keV (Long 
et al. 1990). 

We start with the classical Bethe-Mott-Massey for- 
mula (e.g., Powell 1976), which for the K-shell ionization 
from an ion with i bound electrons reads 



[Tre 4 /(EB^ K )}Z K b hK ln[4c hK E/B hK ], 



(20) 



where S^k is the K-binding energy (taken positive), E{> 
Bi,k) is the incident electron energy, and Zk = 2 is the 
number of K-electrons, and we treat fr^K and c^k as ad- 
justable parameters. It turns out that the existing exper- 
imental data for K-ionization of Ti by electrons in the 
relevant E range of 6 - 50 keV (Long et al. 1990) can be 
well reproduced by the choice of 



1/c 22 ,k w a + (4 - a)exp[d(E/B 22 ,K - 1)] , 



(21) 



with a = 0.55, d = 0.45, and b 22 ,K = 0.445. For sim- 
plicity, we use the same parameter values for the higher 
shells as well as for highly-ionized cases. Aside from the 
trivial replacements of Zk by Zli etc., we use the binding 
energies computed for each ion by the relativistic formal- 
ism, because the rates are sensitive to those thresholds. 
For 44 Ti, their values in keV are, e.g., B^k = 6.6, £?2,k = 
6.2,B 3 , K = 6.1,S 3 ,li = 1.4,B 4 ,k = 5.9,Xli - 1.3, in 
contrast to i?22,K = 5.0, B22,li = 0.56. 

The same parameter values also reproduce the K-shell 
ionization cross-section data for Ni (Long et al. 1990) 
within the experimental uncertainties, while slightly (up 
to 20 %) over-estimating the limited data points for Mn. 
When considered for Fe, this order of the uncertainty of 
the ionization rates influences the radioactivity under con- 
sideration in this work only at a negligible level. 



where rn is the number abundance of 44 Ti +22 ~\ 

The above procedure is quite well suited for dealing 
with the decays of highly-ionized 44 Ti. In fact, our calcu- 
lations based on a relativistic mean field method (Liber- \ rcc f — n c \~J 
mann et al. 1971) show that the overall errors associated i> 



The recombination rates o/ 44 Ti and 56 Fe: The recom- 
bination rates are given by 



_(K) 



9i,LI 



r ( LI ), 



...) , (22) 
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where qt = 1 — pi are the vacancies of the respective shell 
prior to the recombination. We adopt the well-known for- 
mulae (Stobbe 1930) for the recombination cross-sections 
in the non-relativistic approximation for a pure Coulomb 
field. They are expressed in analytic forms in terms of the 
ratios of the incident electron energies (E > 0) to the 
binding energies of respective shells. The relativistic and 
screening effects can roughly be estimated with the use of 
effective charges, which turn out to be largely insignificant 
for the purpose of the present work. 
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